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SUMMARY 


An analytical model of a five degree of freedom magnetic suspension and posi- 
tioning system is presented. The suspended element is a cylinder which is composed 
of permanent magnet material and the magnetic actuators are air core electromagnets 
mounted in a planar array. The analytical model consists of an open loop representa 
tion of the suspension and positioning system with electromagnet currents as inputs 
and displacements and rates in inertial coordinates as outputs. The uncontrolled 
degree of freedom is rotation about the long axis of the suspended cylinder. 


INTRODUCTION 

This paper describes an analytical model of a five degree of freedom magnetic 
suspension and positioning system. The suspended element is a cylinder which is 
composed of permanent magnet material and the magnetic actuators are air core elec- 
tromagnets mounted in a planar array. The uncontrolled degree of freedom is rota- 
tion about the long axis of the suspended cylinder. The electromagnet array is 
mounted horizontally with the suspended element suspended above the array by repul- 
sive forces. The analytical model consists of an open loop representation of the 
suspension and positioning system with electromagnet currents as inputs and displace 
ments and rates in inertial coordinates as outputs. 


EQUATIONS OF MOTION 

This section presents the equations of motion of the suspended element in gen- 
eral terms. The suspended element, or core, is a cylinder which is composed of per- 
manent magnet material. The core is suspended, by repulsive forces, over a planar 
array of five electromagnets mounted in a circular arrangement. Since only five 
degrees of freedom are being controlled, this represents the minimum number of actu- 
ators. The permanent magnet material provides the necessary magnetic field which 
interacts with the fields of the suspension electromagnets to produce suspension and 
positioning forces and torques. A schematic representation of this system which 
shows the coordinate systems and initial alignment is shown in figure 1. A set of 
orthogonal x, y, and z body fixed axes defines the motion of the core with re- 
spect to inertial space. The core axis system is initially aligned with an orthogo- 
nal x, y, and z system fixed in inertial space. A set of orthogonal x^, y^, 
and z b axes, also fixed in inertial space, define the location of the electromagnet 
array with respect to the x, y, and z system. The x and y axes are parallel 
to the x^ and y^ axes respectively, and the z and z^ axes are coincident. 

The centers of the two axis systems are separated by the distance h. 

The development of the core equations of motion follow basically the development 
of the spacecraft equations of motion presented in reference 1. The assumptions used 
for deriving the equations of motion are: (1) the core is a rigid body, (2) the core 
has negligible products of inertia, and (3) there is no motion about the x axis. 
Referring to figure 1, the x axis is an axis of symmetry so that 
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The momentum of the core in core coordinates can be written as 


where 
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A bar over a variable indicates that it is referenced to core coordinates. The total 
torque, •jTk about the center of mass of the core referenced to core coordinates is 


• • 

{t} - {h} » [!]{«} + {fl} x ai]{o}) 


which can be written as 


{t} - [I]{n} + [0][I]{6} 


where [Q] is the skew symmetric cross product matrix defined by 
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can be written 


{*> - ft} + ft J 


where ftj are the control torques on the core produced by the electromagnets and 
■|t V are external disturbance torques. Rearranging terms in equation (6) results in 


{6} - [if 1 <{t} - [6][ij{o}) 
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From equation (1) and the assumption that - 0, the cross product term becomes 
zero and equation (9) reduces to 


{«} - <1/I c ){i} 


( 10 ) 


The angular rates of the core are obtained by integrating equation (10). The core 
Euler rates can be written as 


{?} - [T E ]{n} 
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( 12 ) 


gives the transformation matrix from core body rates to core Euler rates for a 3, 2, 
1 (z, y, x) rotation sequence. 

The forces on the core, in core coordinates, can be expressed as 


{ p } - « c ({°} + («]{«}) 


(13) 


where m is the mass of the core, 
ity. ^fJ can be written 


[ 0 ] 


is defined above, and 



is core veloc- 


0) - { f J + < 14 > 

where ftJL are the control forces on the core produced by the electromagnets and 
4 f | are disturbance forces. Rearranging (13) results in 



(l/m c ){F} - [fl]{v} 


(15) 


The core translational rates are obtained by integrating (15). The core transla- 
tional rates in inertial coordinates, -{ V V , become 


{v> - 


(16) 
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where 
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gives the vector transformation from core coordinates to inertial coordinates. Be- 
cause of the size of the expression, sin has been shortened to s and cos has 
been shortened to c. The displacement of the core center of mass in inertial coor- 
dinates is obtained by integrating equation (16). Since the core is to be actively 
controlled, the assumption is made that the rates will be small and that products of 
rates can be neglected. Under this assumption, the cross product term in equation 
(15), [n]4vk can be dropped and the result is 


{v} - <l/m c ){F} 


(18) 


MAGNETIC TORQUES AND FORCES 

The equations for torques and forces on a magnetic core which are produced by 
air core electromagnets have been presented and discussed in a number of papers on 
wind tunnel magnetic suspension. For example, see references 2-5. This section pre- 
sents the pertinent equations in general terms. 

The torque on a magnetic core, in a given orthogonal x, y, z coordinate 
system, in a nonuniform, but steady, magnetic field can be written as 

{t} - f (d{M) x {b})c!V (19) 


where the integral is taken over the volume of the core and M is the magnetization 
of the core in amp/meter, B is flux density in Tesla, V is volume in cubic me- 
ters, and torque, T, is in Newton-meters. The force on the core can be written as 

{f} - f (d{M}«V){B}dV (20) 

where the force, F, is in Newtons and V is the gradient operator defined as 


v T - [d/dx a/ay a/az] 


( 21 ) 


Since the size of the core is small relative to the size of the electromagnets and 
the gaps involved, a reasonable approximation of the torques and forces which are 
produced can be obtained by making the assumption that the field and gradient 
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components are uniform over the volume of the core. Under this assumption, the 
torque becomes 


{t} - (Vo1)({m} x {b}) 


( 22 ) 


where Vol is the volume of the core and { B ) is calculated at the center of the 
core. Equation (22) can be written as 


{t} - (Vol) [M]{b} 


(23) 


where [M] is the cross product matrix. The force becomes 
{ F } - (Vol) ({m}.V){b} 

Taking the dot product, {») •V, results in the scalar 
{m}.v - {m} t v - (M x a/5x + M y 3/3y + W^/dz) 

Equation (24) then becomes 


(24) 


(25) 


{f} - (Vol) 


M^B^ax + My5B x /3y + M z 3B x /3z' 
\8l y /dx + MyflBy/ay + M z 3B y /3z 
L^aBg/ax + MyaB z /ay + M z aB z /azj 


(26) 


The notation can be simplified considerably by letting represent dB^/dj . By 

factoring M out as a vector, (26) can be written in the form 


( F > - (Vo1)[3B]{m) 


(27) 


where 
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in simplified notation. From Maxwell's equations, in the region of the core 
V x B - 0 which results in 
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( 29 ) 


B - B 
xy yx 

B xz “ B zx 

B yz “ B zy 

Also V*B - 0 which results in 


(30) 

(31) 


B xx + B yy + B zz “ 0 


(32) 


Equations (23) and (27) are written in terms of a given orthogonal coordinate 
system. In the magnetic suspension and positioning system is defined in core 

coordinates and ibV is defined in inertial coordinates. Therefore, in order to 
calculate torque in core coordinates, has to be transformed into the same coor- 

dinate system. This results in 


{i> - <Vol)[M][T m ]{B} 


(33) 


In order to perform the operations involved in the force calculations, {«> is 
transformed into inertial coordinates 



(34) 


and the resulting force transformed back into core coordinates 


- (Vol)[T m ][3B)[T m ] 'H*> 


(35) 


Also, since the core is made of permanent magnet material which is magnetized along 
the x axis, M can be written as 


{m} - Lm- 0 oj 


(36) 


SYSTEM EQUATIONS 

In this section the magnetic torques and forces are combined with the equations 
of motion of the permanent magnet core to produce an open loop model of the magnetic 
suspension and positioning system. Substituting (33) into (10) results in 
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(37) 


{“} - <1/I c )(vo1<[K)[I b ]{b}) + {ij) 

Substituting (35) into (18) results in 

{v} - (l/n, c > (vol([T m )[SB][T m ]- l {M}) + fj (38) 

This gives the equations of motion in core coordinates in terms of and [SB] , 

the gradients of 4bV. Since the core is assumed to be accurately held to a fixed 
point in translation with only small displacements being allowed about the operating 
point, {bV and [SB] at this point are assumed to be essentially constant. This 
means that 4 bV and [SB] can be represented as functions of coil currents only. 

If -ibV and "[SB] are calculated at the operating point, in inertial coordinates, 
for the maximum current in each coil, then 4 B V can be written as 

<39 > 

where I max is the maximum coil current, [Kg] is a 3x5 matrix whose elements 
represent the values of 4 B V produced by a corresponding coil driven by the maximum 
current and 


{i} T - Lij i 2 i 3 i 4 i 5 J 


(40) 


The gradients can be put in the same form by arranging the elements of [SB] as a 
column vector. This results in 


M - "wW 


(41) 


where M is a nine element column vector containing the gradients of ^B^, and 
[K^g] is a 9x5 matrix whose elements represent the values of -[SB^ produced by a 
corresponding coil driven by the maximum current. Each element of Tsb}’, 


for example 


B xx , can be written in the form 


®xx ^/■'•max^L^xx-K*} 


(42) 


where |JScxJ * s a matr i x containing values of B xx produced by a correspond- 

ing coil driven by the maximum current. An analytical model of the magnetic suspen- 
sion and positioning system can now be put into block diagram form as shown in fig- 
ure 2. This model is nonlinear because of the combination of states resulting from 
the coordinate transformations and is of the form 


x - f(x,u) 


(43) 
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where x is given by 


x T -|n n e e v v v xyzl 
l-yzyzxyz 


and the input u is given by 


u T - ll x I 2 I 3 I 4 I 5 J 


(44) 


(45) 


The system can be linearized about a nominal operating point x Q , u Q by performing a 
Taylor series expansion (ref. 6). Ignoring second order terms, the Taylor series 
expansion about x Q , u Q gives 


x + fix - f(x + 5 , u + fi ) 
o o x o u 


(46) 


f(x c , u c ) + df(x, u)/3x 


fix + 3f(x, u)/3u 
x ,u 


5u 


x ,u 
o o 


The equations for ^computing u Q for a given x Q are presented in the appendix. 
Subtracting out x q results in 


fix - Afix + Bfiu 


(47) 


where A - 3f(x, u)/3x| and B - 3f(x, u)/3u| . These equations can be 

I X i VI I X • Vi 

loo loo 

used to calculate feedback gains to stabilize and control the core about the selected 
operating point. Figure 3 illustrates this approach in block diagram form. 


CONCLUDING REMARKS 

This paper has developed an analytical model of a five degree of freedom mag- 
netic suspension and positioning system. The system consists of a suspended element 
which is a cylinder composed of permanent magnet material and magnetic actuators 
which are air core electromagnets mounted in a planar array. The analytical model is 
an open loop representation with displacements and rates in inertial coordinates as 
outputs and electromagnet currents as inputs. The number of electromagnets used in 
the model is the minimum number of five. However, adding electromagnets should be 
straightforward. The model should be useful in analyses and simulations in the de- 
velopment of control systems approaches and in evaluations of overall systems 
performance . 
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APPENDIX 


INITIAL CONDITIONS 

This appendix develops a linearized model of the Magnetic Suspension And Posi- 
tioning System shown in figure 2 by assuming that the suspended core is held fixed at 
a given initial operating point. The linearized model is then used to calculate ini- 
tial conditions on the inputs. The assumptions made in developing the model are that 
the suspended core is held fixed at an initial nominal operating point which can be 
anywhere within 360 degrees in yaw (# z ) but will be zero in pitch (0 y ) . Also, the 
only nonzero disturbance input is assumed to be the weight of the core. Under the 
above assumptions, equation (17) becomes 


t T ml 


-1 


c6 z -s8 z 0 
s6 z c6 z 0 
0 0 1 . 


(Al) 


and 


t T m] - 



s8 z 0 
c6 z 0 
0 1 . 


(A2) 


where $ z is a given fixed yaw angle. From equation (33) the torque in core 
coordinates is 


{t> = (Vol)[M][T m ]{B} 


(A3) 


Substituting (A2) into (A3) and carrying out the matrix operations results in 


(t) - (Vol) 


0 

- M i B z 

_M*(-s* z B x + c* z B y ) 


(A4) 


Factoring out Mjj and substituting from equation (38) for B gives 



(Vol)M-/I 


max 


-L k zJ 

-s0 2 IkJ x + c ^z L^yJ 



(A5) 
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Equation (A5) can be written in the form 


0) - (Vol)M-/I max [KT]{l} 


(A6) 


where [KT] is defined as 


[KT] 


0 

■ L K z J 
L-^zLKxl + o'rl'SrJ. 


(A7) 


Substituting (Al) and (A2) into (35) and performing the matrix operations results in 




(Vol)M- 


c2 ^z B xx + c ^z s ^z B xy + c ^z s ^z B xy + s ^z B yy 
' c ^z s ^z B xx ‘ s tf z B xy + c *z B xy + c ^z s *z B yy 
c *z B xz + s *z B yz 


(A8) 


where the relation B X y - By X has been used. From equation (40), equation (A8) can 
be written as 




(Vol)M-/I 

x max 


c 2 *zK*J + 2c »z s *zKyJ + ‘^zL^yJ 
■ctf z stf z [^K xx j + (c 0 Z - s ^zjL^xyJ + z^^zl^yyl 

“zlAzJ + s# zM 




(A9) 


Equation (A9) can be written in the form 


{f} - <Vol)M-/I max [KF]{l} 


(A10) 


where [KF] is defined as 


[KF] 


c ^zL^xxJ + ^ c ^z s ^zL^xyJ + s ^zl^yyJ 
’ c ^z s ^zL K xxJ + ( c2 ^z ‘ s2 ^z)L K xyJ + c9 z s6 z\?yy\ 
c ^zL^xzJ + s ^z\]^yz\ 


(All) 
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Dropping the zero terra, which is T x , from (A6) and combining with (A10) results in 


S 


(Vol)M-/I 

x max 


[KT]' 


L [KF]J 




(A12) 


where < is a 5x1 column vector and 

IfJ 


’[KT]' 


L(KF]J 


is a 5x5 matrix. 


The solution for the initial currents is then 

•1 

I I K1 I I 

- I /(Vol)M- 
l ) max x 


fi\ 

If) 


(A13) 


‘UKF]J 

As mentioned earlier, the only nonzero element of 


—l is the weight of the core. 
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magnetic suspension & positioning system. 



Figure 2. -Block diagram of analytical model of 
magnetic suspension & positioning 
system. 



Nonlinear System x(t) 

x = f(x,u,t) 
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Figure 3.-Approach for controlling nonlinear system about 
nominal operating point x 0 .u 0 . 
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